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Abstract 

A model for Lai_ x Sr x Mn03 which incorporates the physics of dynamic Jahn- 
Teller and double-exchange effects is presented and solved via a dynamical 
mean field approximation. In an intermediate coupling regime the interplay 
of these two effects is found to reproduce the behavior of the resistivity and 
magnetic transition temperature observed in Lai_ x Sr x Mn03. 
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In this note we present and analyse a model which captures the important physics of 
the "colossal magnetoresistance" materials La!_ x A x Mn03 (here A is a divalent element 
such as Sr or Ca). In the interesting doping range 0.2 < x < 0.4, Lai_ x A x Mn03 is a 
ferromagnetic metal at low temperature, T, and a poorly conducting paramagnet at high T; 
the ferromagnetic-paramagnetic transition occurs at an x dependent transition temperature 
T c (x) ~ 300K and is accompanied by a large drop in the resistivity [1]. The "colossal 
magnetoresistance" which has stimulated the recent interest in these materials is observed 
for temperatures near T c (x) [2]. 

Some aspects of the physics of Lai_ x A x Mn03 are well established [1]. The electronically 
active orbitals are the Mn d-orbitals and the mean number of d electrons per Mn is 4-x. 
The cubic anisotropy and Hund's rule coupling are sufficiently large that 3 electrons go into 
tightly bound d xy , d yz , d xz core states and make up an electrically inert core spin S c of 
magnitude 3/2; the remaining (1-x) electrons go into a band of width ~ 2.5 eV made mostly 
of the outer-shell d x 2_ y 2 and d 3z 2_ r 2 orbitals [Q]. The outer shell electrons are aligned to the 
core states by a Hund's Rule coupling J H which is believed to be large [1]. 

The large value of Jh means that the hopping of an outer shell electron between two 
Mn sites is affected by the relative alignment of the core spins, being maximal when the 
core spins are parallel and minimum when they are antiparallel. This phenomenon, called 
"double exchange" f|], has been widely regarded |||6| as the only significant physics in the 
regime 0.2 < x < 0.5. However, we have previously shown [|7| that double exchange alone 
cannot account for the very large resistivity of the T > T c phase || or for the sharp drop in 
resistivity just below T c . We suggested that the necessary extra physics is a strong electron- 
phonon coupling due at least in part to a Jahn- Teller splitting of the Mn d 4 state in a cubic 
environment. The cubic-tetragonal phase transition observed for < x < 0.2 is known to be 
due to a frozen- in Jahn- Teller distortion with long range order at the wave vector (n, n, n) || . 
We proposed that for x > 0.2 and T > T c (x), slowly fluctuating local Jahn- Teller distortions 
localize the conduction band electrons as polarons. The interesting physics issue is then how 
the polaron effect is "turned off" as T is decreased through T c , permitting the formation of 
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a metallic state. Our picture is as follows. The competition between electron itineracy and 
self-trapping is controlled by the dimensionless ratio of the Jahn- Teller self-trapping energy 
Ej_t and an electron itineracy energy which may be parametrized by an effective hopping 
matrix element t e ff- When Ej_T/t e ff exceeds a critical value we expect a crossover from a 
Fermi liquid to a polaron regime. In models with both double exchange and a large Ej_ T , 
an interesting interplay may occur because t e ff is affected by the degree of magnetic order 
and conversely. As T is increased from zero, the spins begin to disorder. This reduces t e ff 
which increases Ej_t / t e ff so phonon effects become stronger, further localizing the electrons 
and reducing t e ff and thereby the effective ferromagnetic coupling. 

To investigate this quantitatively we consider the model Hamiltonian H eff = H e i + Hj_ T 
with 

Hei = - E ty'dLdjb. + JnEt 4^cU + h • S c /S c (1) 

ija i,a,a 

and 

Hj-t = gE d LQ ab (j)<V + kE Q 2 0)- (2) 

jao- j 

Here d^ CT (i) creates an outer-shell d-electron of spin a in the a orbital on site i. The local 
lattice distortions which cause the Jahn- Teller splitting transform as a two-fold degenerate 
representation of the cubic group which we parametrize by a magnitude r and an angle (f). 
They couple to the electron as a traceless symmetric matrix Q = r(cos(0)r z + sin(0)r x ). The 
electron-phonon coupling is g and the phonon stiffness is k. The external magnetic field is 
h; for simplicity we have coupled it to the core spin only. In the phonon part of H eff we have 
neglected intersite terms and also cubic and higher nonlinearities. In the electronic part 
of H cff we have neglected on-site Coulomb interaction effects; these will be important for 
higher energy properties of spectral functions but will only affect the low-energy properties 
we consider here by renormalizing parameters such as ty b . 

To solve H eff we introduce further simplifications. We take Jh — > oo. Because we are 
interested in phenomena at temperatures of order room temperature, we assume the phonons 
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and the core spins are classical. We allow for magnetic order but assume that there is no 
long range order in the lattice degrees of freedom. To solve the electronic problem we use 
the "dynamical mean field" approximation which becomes exact in a limit in which the 



spatial dimensionality d — > oo JT0|. Then, the free energy may be expressed in terms of a 



space-independent "effective field" G e ff(u;) via 

Z = J rdrd0dfiexp[-tr 2 /2T + TrlnftGetr" 1 + Ar • f + J U S C ■ a] +h • fi] (3) 

Here Q is the direction of S c and t = D/4 (D is the full bandwidth, so from || one estimates 
t w 0.6eV). The dimensionless electron-phonon coupling constant A = g/Vkt. G e ff(w) 
is a tensor with orbital and spin indices; it obeys a self-consistency condition whose form 



depends upon the lattice whose d — > oo limit is taken [10]. We have used the Bethe lattice 



equation, which corresponds to an underlying band structure with a semicircular density of 



states with D = 2Trt . The self consistent equation is [pLO | 



G ef f _1 H =w-/i-Tr[tGt]/2 (4) 

where G = 91nZ/5G e ff _ . Because we assume there is no long range order in the lattice 
degrees of freedom, we take G e ff to be the unit matrix in orbital space. We have used 
two methods for treating the spin part of the problem. In the direct integration method, 
one solves Eq. (4) by performing the integrals over the angle and phonon coordinates 
numerically. In the projection method, one quantizes the electron spin on site i along an 
axis parallel to SJ. and retains only the component parallel to S*. The Jh term then drops 
out of the Hamiltonian but as shown previously 0,0], one must multiply ty by the double 



exchange factor qy = cos(#ij/2) = y (1 + Sj. ■ S J c )/2. Within mean field theory one may 
replace q^ by q = (l + m 2 /2)/v^2, where m = (SJ.)/S C is determined self consistently via m = 
-T<9/dh[coth/3(Jm + h) - (/3(Jm + h))" 1 ) and as shown previously 0, J = (l/2y/^)dbaZ/dt 



with Z evaluated at q = y (1 + m 2 )/2. The resulting d = oo equations involve a G c fr which 
is a scalar and a numerical integral over the phonon coordinate only. Because t enters the 
mean field equations only as an energy scale, it is necessary only to solve the resulting mean 



field equations once at each T and A to yield a Z(T/t, A 2 /t); the q dependence and hence the 
magnetic properties may be found by scaling t — > qt. The two approaches give very similar 
results for the magnetic phase boundary and the phonon contribution to the resistivity, but 
the direct integration approach gives also the spin disorder contribution to the resistivity. 
The main difference between the two calculations is that projection method leads to a first 
order magnetic phase transition for A ~ 1. 

We now discuss the solutions. Several soluble limits exist. At A = there is a second 
order ferromagnetic transition at a T c (x) which is maximal for the half filled band (T c (0) = 
0.1 7t) and decreases as the band filling is decreased. For T > T c there is spin disorder 
scattering which (in the mean field approximation used here) is temperature independent 
and of small magnitude. This scattering decreases below T c as discussed previously ]5] 

In the limit T — > 0, ground state is a fully polarized ferromagnet for all A. The phonon 
probability distribution P(r) = J d</>dfie~ tr2 / 2T+Trln [ tGeff 1+Ar ' f l is sharply peaked about the 
most probable value r = r*. For A < A c (x), r* = and the ground state is a conducting Fermi 
liquid. For \ cS > A c , r* > 0, implying a frozen-in lattice distortion and, if r* is large enough, 
a gap in the electronic spectrum. For x = 0, A c = 1.08... and the transition is second order, 



with r* linear in ^A — A c but very rapidly growing, reaching the point r* = 1 at A = 1.15. 
For r* > 1, r* becomes linear in A and a gap appears in the spectral function. For x > 
the transition is first order, involves a jump to a nonzero r*, and occurs at a A c > 1.08. A 



detailed discussion of the spectral functions and transitions will be presented elsewhere fl2 
The increase of A c with x is due to the increased kinetic energy per electron. Note that 
the double exchange effect means that the kinetic energy is maximal in the fully polarized 
ferromagnetic state. For uncorrelated spins, the kinetic energy is smaller by a factor of \[2 
and A c smaller by a factor of 2 1//4 . In other words, there is a regime of parameters in which 
the electron-phonon interaction is insufficient to localize the electrons at T = but sufficient 
to localize them at T > T c (x). 

Another analytically solvable limit is A 3> 1. In this limit an expansion in 1/A may be 
constructed for arbitrary 1/AT, and in the leading order one may evaluate the r integral by 



steepest descents. Here we find a second order phase transition at T c = t/12A 2 separating 
two insulating phases with slightly different gaps. 

We turn now to numerical results, limiting ourselves to x=0 for simplicity. 

The Jh — > oo limit means that the d-bands arising from the outer-shell orbitals are half 
filled, so the chemical potential \x = 0. Eq (4) is solved on the Matsubara axis by direct 
iteration starting with the A = solution. From this solution Z is constructed and (m) and 
(r) are computed. The conductivity is calculated following fllQf ; the requisite G e ff on the 



real axis is obtained by solving Eq. 4 for real frequencies, using the previously obtained 
Matsubara solution to define Z. 

Fig 1. shows the phase diagram in the T — A plane. The solid line is a second order 
transition separating ferromagnetic (F) and paramagnetic (P) regions obtained via the direct 
solution method. The light dashed lines separate regions of weak electron-phonon coupling 
in which dp/dT > from regions of strong electron-phonon coupling in which dp/dT < 0. 
We identify these regions as metal (M) and insulator (I) respectively. The T dependence of 
the dp/dT line below T c is due mostly to the temperature dependence of the magnetization. 
Increasing A decreases T c ; the variation is particularly rapid in the crossover region A ~ 1, 
consistent with the very rapid dependence of r* on A mentyioned above. Here also the 
magnetic transition calculated via direct integration becomes more nearly first order. The 
projection method leads to a region of two-phase coexistence for 0.92 < A < 1.1. This is 
shown on fig 1 as the area between the heavy dotted line and the solid line. The different 
behavior of the two models suggests that in the crossover region the order of the transition 
is sensitive to the approximations of the model. Other physics, not included here, will also 
tend to drive the transition first order. We mention in particular anharmonicity in the 
elastic theory, which will couple the Jahn- Teller distortions to the uniform strain, and also 
the conduction electron contribution to the binding energy of the crystal, which produces 



the observed volume change at T c JTjj] . 

The inset to fig 1 shows the average of the square of the lattice displacement. This 
would be measurable in a scattering experiment sensitive to rms oxygen displacements. In 
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the classical model used here, r 2 — > r 2 + T as T — > 0. One sees that for intermediate 
couplings the high temperature state has a non-zero extrapolation to T = while the low 
T state has a vanishing extrapolation, while for larger couplings both sides of the transition 
have non-zero but different extrapolations. 

Fig 2 shows the temperature dependence of the calculated resistivity for several different 
values of A. At small A and T > T c , p is small and has a T- independent piece due to the 
spin disorder and a T-linear piece (difficult to percieve on the logarithmic scale used in fig 
2), due to electron-phonon scattering. As T is decreased through T c , p drops as the spin 
scattering is frozen out and the phonon contribution changes slightly. For larger A a gap 
opens in the electron spectral function at T > T c and the resistivity rises as T is lowered 
to T c and then drops sharply below T c , as the gap closes and metallic behavior is restored. 
Finally, at still stronger coupling, insulating behavior occurs on both sides of the transition, 
although there is still a pronounced drop in p at T c . 

Fig 3 shows the magnetic field dependence of p for A = 1.12, demonstrating that in 
this region of the phase diagram the "colossal magnetoresistance" phenomenon occurs. The 
magnetic field scale is too large relative to experiment (as is the calculated T c ), but is very 
small in comparison to the microscopic scales of the theory. 

The series of resistivity curves presented in fig 2 bears a striking resemblance to measured 
resistivities on the series Lai_ x A x Mn03. We have already noted that the calculation, which 
neglects long range order, gives the generic behavior at any carrier concentration. We identify 
the experimental doping x with the relative strength of the electron phonon interaction, A 
because increasing x increases the kinetic energy per electron. With this identification the 
results are consistent with the observed variation of T c and p with x, and also with the 



opening of a gap observed [|T3[ in the optical conductivity. Note also that as T moves 
through T c , the effective ferromagnetic J should drop by about 10 percent, producing a 
perhaps observable shift in the position of the zone boundary magnon. 

The present theory is consistent with a recent study of a La(Pr, Y)Ca.3Mn03 series 
of compounds 0]. The substitution of Pr,Y for La decreases the effective d-d overlap, 



decreasing t and increasing A. Experimentally, it results in a shift of T c to lower temperatures 
and an increasing resistivity anomaly, as found in the calculation. The observed first order 
transition also occurs in one of the mean field theories we have considered. 

In summary, we have presented a solution of a model describing the Jahn- Teller and 
double exchange physics, and have shown that it accounts naturally for the existence of a 
high-T insulating phase, the dramatic changes of resistivity at T c , and the extreme sensitivity 
to magnetic field. In obtaining this behavior the interplay of polaron and double exchange 
physics is essential. We have solved the d = oo equation appropriate to Eq. 3 with Jh = 
and found metallic behavior (with p and r 2 linear in T) at temperatures greater than the 
T = gap and insulating behavior at temperatures less than the T = gap. A more detailed 
exploration of the MF theory allowing uniform or staggered ordering of the lattice distortions 
and varying carrier concentration will be needed to study the structural transition at low 



x and the "charge ordered" phase flTT| , |T5|l at x « .5. These calculations however must also 
include the intersite phonon coupling. 

We acknowledge stimulating discussions with G.Aeppli, S-W. Cheong, A. Georges, H. 
Hwang, B. G. Kotliar, H. Monien, A. Ramirez, T. M. Rice, M. Rozenberg, P. Schiffer and 
R. Walstedt. We are particularly grateful to P. B. Littlewood, who stimulated our interest 
in the problem, collaborated in the early stages of this work, and has been a continuing 
source of help and encouragement. A. J. M. thanks the Institute Giamarchi-Garnier and 
the Aspen Center for Physics and B. I. S. the Ecole Normal Superiure for hospitality. R. M. 
was supported in part by the Studienstiftung des Deutschen Volkes. 
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Figure Captions 

Fig. 1: Phase diagram. Solid line: ferromagnetic T c as a function of electron-phonon 
coupling, A, calculated by direct integration method. The area enclosed by the solid line 
and the heavy dashed line is the region of metastability found in one formulation of mean 
field theory. Light dotted lines: metal-insulator crossover obtained from calculated re- 
sistivities. Regions labelled as PM (paramagnetic metal), FM (ferromagnetic metal), PI 
(paramgnetic insulator) and FI (ferromagnetic metal) according to the value of the magne- 
tization and dp/dT. Inset: square of average lattice distortion plotted vs temperature for 
A = 0.71 (lowest), 0.9, 1.05, 1.12, 1.2 

Fig. 2: Resistivity calculated by direct integration method plotted versus temperature 
for different couplings A. Heavy solid curve (top) A = 1.2, heavy dotted curve, A = 1.12, 
heavy dot-dash curve A = 1.05, light solid curve A = 0.95, light dashed curve A = 0.85, light 
dot-dashed curve (bottom) A = 0.71. 

Fig 3. Magnetic field dependence of resistivity calculated by direct integration method 
for A = 1.12 and magnetic field h as shown. Note h=0.01t corresponds to 15 Tesla if 
t = 0.6eV and S c = 3/2. 
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